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Abstract 

We study light-like polygonal Wilson loops in three-dimensional Chern-Simons 
and ABJM theory to two-loop order. For both theories we demonstrate that the 
one-loop contribution to these correlators cancels. For pure Chern-Simons, we find 
that specific UV divergences arise from diagrams involving two cusps, implying the 
loss of finiteness and topological invariance at two-loop order. Studying those UV 
divergences we derive anomalous conformal Ward identities for n-cusped Wilson 
loops which restrict the finite part of the latter to conformally invariant functions. 
We also compute the four-cusp Wilson loop in ABJM theory to two-loop order 
and find that the result is remarkably similar to that of the corresponding Wil- 
son loop in A/" = 4 SYM. Finally, we speculate about the existence of a Wilson 
loop/scattering amplitude relation in ABJM theory. 
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1 Introduction and conclusions 



Wilson loops are the central non-local observables in any gauge theory and thus of intrinsic in- 
terest. In 3d Chern-Simons theory they are the principal observables and topologically invariant 
with exactly known correlation functions [Ij in the Euclidean (or Wick rotated) theory. This 
exact result is analytic in the inverse Chern-Simons parameter k and perturbative studies in a 
loop-expansion of the effective coupling constant 1/k can reproduce the exact topological and 
finite result to the first orders [2,3J, modulo regularisation subtleties leading to or not leading 
to an integer shift of k (for a review see [4J). Wilson loops in Minkowski-space with cusps and 
light-like segments, however, display particularly strong divergences in 4d gauge theories and 
seem to not have been considered in the 3d Chern-Simons literature before. 

In this paper we study such light-like Wilson loops with cusps of polygonal shape in pertur- 
bation theory up to the next-to-leading order in 1/k. We do this for both the pure 3d Chern- 
Simons theory as well as its conformal J\f = 6 supersymmetric extension known as Aharony, 
Bergman, Jafferis, Maldacena (ABJM) theory [5J. The A/" = 6 ABJM theory is built upon an 
SU{N) X SU{N) gauge symmetry which allows for a planar -^^ oo limit with A = N/k held 
fixed, where k is the common absolute value of the Chern-Simons parameters of the two SU {N) 
subgroups. In this limit the ABJM theory is conjectured to be dual to type IIA string theory 
on AdS4^ X CPs, representing an exact gauge-string duality pair very similar in nature to the 
well studied 4d A/" = 4 super Yang Mills/ AdS^ x string duality pair. Supersymmetric Wilson 
loops in ABJM theory have been defined in [6-8J for the 1/6 BPS and recently in [9J for the 1/2 
BPS case. The correlators for these loops of Euclidean, circular geometry are moreover known 
exactly in terms of a supermatrix model [10] using localisation techniques [llj. Our motivation 
to consider polygonal light-like Wilson loops in the 3d Chern-Simons ABJM gauge theory stems 
from the Wilson loop/scattering amplitude duality in A/" = 4 super Yang-Mills. This duality was 
discovered in the dual AdS^ x string picture at strong gauge coupling in [12j and shown to 
exist also in the weak coupling regime [13-15J with profound consequences on the symmetries 
of these correlators leading to a dual super conformal [16j respectively Yangian symmetry [17j of 
scattering amplitudes, for reviews see [18,19j. Moreover, there are many structural similarities 
of the 3d A/" = 6 superconformal ABJM theory to N' = 4 super Yang-Mills, most notably the 
emergence of hidden integrability [20-23J, (for reviews see [24-31J) in the planar limit for the 
spectral problem of determining anomalous scaling dimensions of local operators [32,33j. 

Given these insights the question arises whether there could also be such a scattering ampli- 
tude/Wilson loop duality in the ABJM theory. Scattering amplitudes in the ABJM theory have 
been analysed by Agarwal, Beisert and McLoughlin [34J who in fact studied more general mass 
deformed superconformal Chern-Simons theories with extended supersymmetries at the one-loop 
order. There a vanishing result for the four-point one-loop amplitudes in the ABJM theory was 
found and the authors speculated whether the two-loop scattering amplitudes in A/" = 6 Chern 
Simons (ABJM theory) could be simply related to the one-loop M' = 4 Yang-Mills amplitudes. 
The main result of our paper is that this picture is consistent at least up to the two-loop order: 
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We observe a cancellation of the one-loop graphs for null polygonal loops and find that the four- 
cusp Wilson loop at the two-loop order is of the same functional form as the four-point MHV 
amplitude of A/" = 4 up to constant numerical terms. 

Specifically we calculate the expectation value of the n-cusped Wilson loop operator 

{Wn):=j^{0\TTVexp(^t^^A^dz>'^ |0) (1.1) 

in the planar limit ^ for light-like polygonal contours C in pure Chern-Simons and ABJM theory, 
see Appendix A for conventions of the Lagrangian and the path ordering. 

The contour of the n-sided polygon C is given by n points {i = 1, ...,n) where we parametrise 
each edge Cj via 

Ziisi) = x'^+p'^Si, = - xf (1.2) 

where Si G [0, 1 ] and C = CiU...U C„. The calculations are performed in 3-dimensional Minkowski 
space with metric 

r/^,, = diag(l, -1, -1) , (1.3) 

and the segments of the contour are light-like, i.e. pf = 0. Note that due to the light-like contour 
there is no difference in ABJM theory between the standard loop operator (1.1) above and the 
1/6 BPS supersymmetric loop operator of [6J, as the terms in the exponential coupling to the 
scalar s drop out. 

Wilson loops in Chern-Simons theory are usually defined with a framing procedure [1-3J which 
may be thought of as a widening of the Wilson line to a ribbon. This is necessary in order to 
define an integer twisting number of the individual loop and acts as a particular point-splitting 
regulator for collapsing gauge field propagators in perturbation theory while preserving the topo- 
logical structure of the theory. Here we refrain from framing our loops as we do not encounter 
the problem of collapsing gauge field propagators due to the piece-wise linear structure of our 
loops. Moreover, the ABJM theory is not topological due to metric dependent interactions in 
the matter sector, so that there is no need for framing from that perspective either. Instead 
we regulate our correlators by the method of dimensional reduction which has been tested to 
the three-loop order in pure 3d Chern-Simons to yield a vanishing /3-function and to satisfy the 
Slavnov- Taylor identities [35j. Here the tensor algebra is performed in 3 dimensions to obtain 
scalar integrands and then the dimension of the integrations are analytically continued. 

The outcome of our computations at one-loop order in pure Chern-Simons and ABJM theory 
is that as claimed 

(W^4)i-ioop = (analytically) , (W^6)i-ioop = (numerically) . (1.4) 
^I.e. we take the limit N,k ^ oo, ^ finite. 
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Moreover, conformal Ward identities force (Wn) i-ioop to depend only on conformally invariant 
cross ratios of the (xj — and for n = 4 and 6 we show that these functions vanish. 

At the two-loop order we computed the tetragonal Wilson-loop in pure Chern-Simons 
and TV = 6 superconformal Chern-Simons (ABJM theory). The result in dimensional reduction 
regular isat ion with c? = 3 — 2e for the correlator in pure planar Chern-Simons reads 



cs 



ln(2)- 



-^?3/i')'^. ^ 1^(2)^^!^ + const. 



2e 



2e 



+ o[{ 



1.5) 



We remark that this result displays a breakdown of finiteness and topological nature of the 
light-like four-cusp Wilson loop in 3d Chern-Simons at the two-loop order due to divergences^ 
associated to two cusps at a light-like distance, see section 3.2, For the same correlator in the 
ABJM theory we find 



ABJM 



(2e) 



+ 



+ (1^1+ const. 



;i.6) 



This is indeed of the same functional form as the one-loop result in A/" = 4 super Yang-Mills 
theory, where one has [13j 



Ar=4 SYM 



X 



13 



24 



+ const. 



;i-7) 



It would certainly be interesting to explore this relationship also beyond four cusps. Also a 
two-loop computation of the four-particle scattering amplitudes in ABJM theory would be very 
desirable in order to compare to our result. 



From the string perspective the scattering amplitude/Wilson loop duality in the AdS^/CFT^ 
system arises from a combination of bosonic and fermionic T-dualities under which the free 
AdS's X superstring is self-dual [36,37J. Hence, for the existence of an analogue duality in 
ABJM theory one would require a similar self-duality of the AdS^ x CP3 superstring under a suit- 
able combination of T-dualities. This problem was analysed by Adam, Dekal and Oz in [38j with 
a negative outcome: The Green-Schwarz a-model is not self-dual under bosonic T-dualtities 
in the transverse AdS^ directions combined with fermionic ones within the framework of the 
fermionic Buscher dualisation procedure employed in [36j. However, the analysis of [38j started 
from a partially K-symmetry gauge fixed formulation of the AdS^ x CP3 superstring in terms 
of a supercoset cr-model [39,40j. This restriction was overcome in [41j where, building upon a 
complete superspace formulation of AdS^ x CP3 [42], again the non-existence of a T-self-duality 



^Regularisation is known to be a highly subtle issue in 3d Chern-Simons theory. In particular one could 
compute this correlator using a framing procedure adapted to light-like curves, which would be very interesting. 
Wc thank N. Drukker for discussions on this point. 
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of the AdS^ x CP3 superstring was found. Interestingly however, a very recent paper [43j has 
uncovered a Yangian symmetry of tree-level amplitudes in ABJM theory pointing towards inte- 
grability of the latter. Moreover, the authors argue about the possibility of a self-duality of the 
AdS^ X CP3 superstring upon T-dualizing also along the CP3 directions, which could provide a 
loop-hole for a scattering-amplitude/light-like Wilson loop duality of the ABJM theory. In order 
to settle this question a two-loop scattering amplitude computation in ABJM theory would be 
very desirable in order to compare to our result ( 1.6) . 

Our paper is organised as follows: In section 2 we perform the one-loop computation in pure 
Chern-Simons and ABJM theory for the tetragon and the hexagon Wilson loop. Section 3 is 
devoted to the two-loop problem of the tetragon in pure Chern-Simons theory where we compute 
all relevant graphs. Then, in section 4, we perform an independent check of our results by 
deriving anomalous conformal Ward identities for the Wilson loop. The latter can be generalised 
to an arbitrary number of points. In section 5 we include the matter diagrams arising in ABJM 
theory at the two-loop order and combine the results for the tetragon to obtain our final result 
(1.6). In the appendix we collect a detailed account of our conventions and give the technical 
details of the computation of the two-loop graphs using Mellin-Barnes techniques. 

Note added 

After publication of this paper the works [44J, [45j appeared which report on a two loop calculation 
of four-point scattering amplitudes in ABJM theory. Most interestingly, the results coincide 
precisely with the divergent and finite pieces of our Wilson loop computation in (1.6) up to a 
constant. In the published version of this article there was an erroneous sign in (5.11) leading to 
a spurious sign difference between the Wilson loop and the scattering amplitude. We thank the 
authors of [45j for pointing this out to us. 

2 One loop: Chern-Simons and ABJM theory 

In this section we consider the one-loop expectation value of polygonal Wilson loops with n 
cusps. We would like to consider kinematical configurations for which all non-zero distances 
satisfy —xfj > 0, such that the result for the Wilson loops will be real (In particular, this allows 
us to drop the ie prescription of the propagators). For n odd, however, it is impossible to find 
vectors that lead to such configurations. For this reason, we will only discuss n even. 

At one-loop level, we only need terms quadratic in the expansion of the Wilson loop operator, 
and the free part of the action. Therefore, at one loop order, the expectation value of (1.1) in 
ABJM theory coincides with the one in pure Chern-Simons theory. 

The expectation value at one loop is a sum over all possible diagrams where the propagator 
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stretches between edges i and j, 

where the domain of integration is given by dsi dsj for i 7^ j, dsi J^' dsj for i = j and 
z{si) = dz{si)/dsi = Pi and where we have introduced a normahsation factor for later convenience. 
Here and throughout the paper, we absorb the dimensional regularisation scale [fi'^Y into k and 
only display it explicitly in our final results. Using the Chern-Simons propagator in the Landau 
gauged 

U.)mn (^) (y)) = ^n^lSnkl ( ^ ) e,., ^^^^ , (2.2) 



and plugging in the expressions (1.2) for 2j, we obtain 



e{pi,Pj,PiSi - PjSj + Y!k=jPk 



lij = I dsidsj — J, (2.3) 



2 



•^ij^i^J -^1+1,3^1^3 -^ij+l^i^J -^i+l^i+l^^^J ) 

where xj^^ = (xj — XjY, Sj = 1 — Sj and e(a, 6, c) = e^ka^Vc^ . We can immediately see that in this 
gauge li^i and h^i+i vanish due to the antisymmetry of the e tensor. This corresponds to diagrams 
where the propagator ends on the same edge or on adjacent edges, as shown in Figures l[b) and 
1(c), respectively. Therefore we only need to keep diagrams of the type shown in Figure 1(a). 
The latter are manifestly finite in three dimensions and therefore we set (i = 3 in the remainder 
of this section. 

2.1 Tetragon 

As explained above, in the Landau gauge, the only non-vanishing contributions to (2.1) for the 
tetragon are /31 and 1^2- Setting li = 3, they are given by 



'31 



-e(Pi,P2,P3) / dsids-i- 2 - - ^ 2 (2-4) 



^We drop the ie prescription in the propagator, since we consider kinematical configurations with —xf, > 0. 
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Figure 2: One-loop contributions Ii+2,i and /i+3,i to the hexagonal Wilson loop. 



and 



^42 = -e(P2,P3,P4) j dS2dSi-^ 



-X 



"i 7. 7^ 

24 



S2S4 - xl^S2SiY/'^ 



(2.5) 



Taking into account that we have a closed contour, i.e. YliiPi — 0' "^^ '^^'^ write 6(^2,^3,^4) = 
— e(p2,P3,Pi) = ~e(Pi,P2,P3) and thus the contributions from the two diagrams cancel each other 



+ 



/3I + /42 = . 



(2.6) 



We will see in section 4 that this result is compatible with the restrictions imposed by conformal 
symmetry. 



2.2 Hexagon and higher polygons 

For the hexagon there are two different non- vanishing types of contributions, /i+2,i and h+zAi as 
shown in Figure 2, The former appears in six orientations, i = 1 ... 6 (with the convention that 
i + 6 = 2), while the latter appears in three orientations, i = 1, 2, 3. 
Specialising the general formula (2.3) to these cases we have 

T f\ , ^{Pi+2,Pi,Pi+l) „^ 

h+2,i = / dsi+2dsi _ _ 2 = 2 2 W9 y^-'i 



/ — — 2 — 2 2 ^3/2 



and 



r _ f\ r eiPi+3,Pi,Pi+i+Pi+2) „^ 

-'i+3,i — / O'Si+30'Si- 2 - 2 - 2 2 \3/2 ' V^'"/ 

We checked numerically for various non-symmetric hexagon configurations that the sum over 
all diagrams vanishes, 

6 

(W^6)^'^cx5^/,, = 0. (2.9) 

i>j 

Although we do not yet have an analytical proof for generic kinematical configurations, we can 
show that (2.9) is true for special configurations, as we will see presently. 
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(a) 



(b) 



I 

(c) 




(d) (e) (f) 

Figure 3: Planar two-loop topologies appearing in the polygonal Wilson loop in CS theory. 
Diagrams where one propagator is connected to a single edge or to two adjacent edges vanish in 
our gauge and are not displayed. 

Consider the configuration where opposite edges are anti-parallel, i.e. Pi = —pi+s- From (2.8) 
we see that /i,i+3 = due to the antisymmetry of the e tensor. Furthermore, taking into account 
that for this configuration we have x^j_^2 = ^?+3,i+5; easy to see from equation (2.7) that 
the integrands of /i,i+2 and Jj+s^j+s are the same. Finally, using ^jPi = one can see that the 
Levi-Civita symbols produce a differing sign, such that 




(2.10) 



i.e. the contributions coming from those diagrams cancel pairwise, and we arrive at equation 
(2.9), in the specific anti-parallel kinematical configuration pi = —pi+3. 

It is tempting to speculate that all n-cusped Wilson loops vanish at one-loop order in Chern- 
Simons theory. 

3 Two loops: Chern-Simons theory 

In the this section we calculate the two-loop contributions to the tetragonal light-like Wilson loop 
in pure Chern Simons theory. The results are consistent with the anomalous conformal Ward 
identity to be discussed in section 4, 

Expanding the Wilson loop to quartic order, see (A. 12), and performing Wick contractions 
leads to the topologies shown in Figure 3, We are taking the planar limit and therefore drop all 
non-planar graphs. Moreover, all diagrams where one propagator is connected to a single edge 
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or adjacent edges vanish in our gauge for the same reason as at the one- loop order and are not 
displayed. 



3.1 Ladder diagrams 

Let us begin by computing diagrams of ladder topology as shown in Figure 3d, There are two 
different orientations of this diagram, and it is easy to see that they give the same contribution. 
Taking into account this factor of 2, we have that the contribution of the ladder diagrams is 



)■(»)■■• 



Wfler = 2(f) /ladder(x?3,X^4), (3-1) 



where 



J ( 2 ™2 N _ I J ^(^ii Zu Zj Zj) e{Zj, Zk, Zj Zk) , . 

-'ladderl.-^lS) -^24/ ~ / '^^i,j,k,l ^ ^ d , , ,„-,d • I'-'-^J 



\-\Zi-zir\^ [-[Zj-Zk) 



2 



The integral is finite and may be calculated for d = 3 

/ladder(x?3,^L) = ^ / dsj dsJ dsj ds,- ^13^24(^13 + ^24) + 0{e) . 

^ Jo Jo Jo Jo [xi^SiSi + X24SiSi\ 2 [xfgSjSfc + X24SjSk\ 2 

(3.3) 

We computed this integral by first carrying out some of the parameter integrals and then deriving 
a differential equation for it, which could be solved. The result is remarkably simple. 



-^ladder(3;i3) 2^24) — ^ 



X 1 



InM^l+vr^ 

24, 



0(e) . (3.4) 



Including the prefactors and dropping 0(e) terms, the contribution to the Wilson loop is 



3.2 Vertex diagrams 



(2) 



ladder 



A^\^ 1 



X24 



(3.5) 



The diagrams with one three-gluon vertex shown in Figures 3al , |3b] and [3c| are obtained by 
contracting the cubic term in the expansion of the Wilson loop in ([A.12p with the interaction 
term of the Lagrangian, 



TC 2 



i>j>k 
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Figure 4: The divergence in the vertex diagram arises from the integration region where si — )■ 
1, S3 — > (and (3i — > 0,/33 — 0), see equation (3.9). 



where 



- [ dz^dz^dzy^^e,^.e.p,e^r [ d'w ^"^ "^^.^j^ '^H'' '"l^ , (3.7) 



and \zi\ = {—zf)^. Here the indices of lijk refer to the edges of the Wilson loop that the 
propagators attach to. The expression can be shown to be antisymmetric under the exchange of 
any two indices, and therefore the only non-vanishing contributions are the ones for i ^ j ^ k. 
As a consequence, topologies 3a and 3b can be discarded. 

Specialising to the tetragon, we have four contributions which are symmetric under xf^ ^ 
and thus it is sufficient to compute one of them 



'321 



,d / \ e(P2,P3,w)e(p2,Pi,w) 

d w / dsi,2,3-r-n71 iJi u (3-8 

\W\'^\W — Zi2\\W — Z32\ 



^^"^ d's,,2,3d'(3,,2,3 (/3l/32/33)'^ 



^d 



2^^^r7-^/3i/33SiS3(xi3 + xl^] 



where the second and third line is obtained by introducing Feynman parameters in the standard 
way and integrating over w. More details may be found in Appendix B.l, Z\ is given by 

^ = -fiihzl^ - /32/334 - /3i/33^?3 

= -Xi3/3iSi (/33S3 + I32S2) - 0:24/3353 {1^282 + /3isi) (3.9) 

One might naively think that this diagram should give a finite answer due to the antisymmetry 
of the e tensors. The result would indeed be finite in the case of smooth contours [2J or contours 
with a single cusp. However, the presence of two cusps gives rise to a region in the integration 
space of Feynman parameters where the first summand in the third line of (3.8) induces a 
divergent contribution. The relevant region of Feynman parameters is si — )• 1, S3 — )■ (and 
/3i — )■ 0, /33 — )■ 0), see equation (3.9), and is illustrated in Figure 4. Due to the presence of three 
independent vectors Pi,P2 and P3 the e tensors do not suppress this region. We find that this 
term produces a 1/e pole in dimensional reduction. The second summand in the third line of 
(3.8) is finite. 
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We separated the divergent and finite pieces using Mellin-Barnes techniques. The details 
can be found in Appendix B.l. We have not computed the coefficients of the and e° terms 
analytically, but we have good numerical evidence that they give the following result: 



'321 



1) 

8 rilY 



21n(2) 



-X 



2 \2e 
24 J 



In^ 



X 



13 



24 



+ 06 + 0(e) 



(3.10) 



where = 8.354242685 ±2-10 ^, see (IB.27p . Taking into account all prefactors and restoring 
the regularisation scale, k — )■ fi~'^'^k, we can write the result, up to terms of order e, as 



(^4)tertex 



4 



(3.11) 



3.3 Gauge field and ghost loops 

It is well known [35j that in the dimensional reduction (DRED) scheme the gauge field loop 
diagrams shown in Figure 3e exactly cancel against the ghost loop diagrams shown in Figure 3f : 



gluon loop ghost loop * 

Details of this cancellation can be found in appendix B.2,, 



(3.12) 



3.4 Result for the two-loop tetragon in CS theory 

Summing up the results (3.5), (3.11) and (3.12) for the tetragon, interestingly, the ln^(x^3/x24) 
terms in the two-gauge-field diagram and the vertex diagram exactly cancel and we obtain 



(2) = _M1 ^ 



1 

4 



ln(2) 



-^1+2 P'r 



+ a6-81n(2) - vr^ 



i=l 



(3.13) 



where fl"^ = ^^T^e""^ and we recall that = 8.354242685 ± 2 ■ 10~^. As we will see in the next 
section, the cancellation observed here that led to the finite part of (3.13) being a constant is in 
fact a consequence of the (broken) conformal symmetry of the Wilson loops under consideration. 



4 Anomalous conformal Ward identities 

The structure of the above results can be understood from conformal symmetry, by deriving 
anomalous conformal Ward identities for the Wilson loops. Here we follow very closely reference 
[46] . 

We would like to use the specific properties of Wilson loops with light-like contours C under 
conformal transformations. The key point is that such contours are stable under conformal 
transformations, i.e. the deformed contour C is also made of n light-like segments. This can be 
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seen as follows. The cusp points Xi form a contour with light-like edges, i.e. xf^j^^ = 0. It is 
obvious that the light-likeness conditions are preserved by translations, rotations, and dilatations. 
Special conformal transformations are equivalent to an inversion — )■ x^/x^ followed by a 
translation and another inversion. Thus it remains to investigate the transformation under 
inversions. Since under the latter xf^ — )■ it is clear that the light-likeness of the 

contour is preserved by all conformal transformations. 

If the Wilson loop (Wn) were well defined in (i = 3 dimensional Minkowski space it would 
enjoy the conformal invariance of the underlying gauge theory and we would conclude that 
{W{C)) = {W {€')). This is indeed the case at one loop order, see section 2. However, as we have 
seen in section 3, starting from two loops, divergences force us to introduce a regular isat ion and 
calculate in d = 3 — 2e dimensions, thereby breaking the conformal invariance of the action. The 
latter leads to an anomalous term in the conformal Ward identities for the Wilson loops, as we 
will see presently. 

The expectation value of the Wilson loop can be written as a functional integral 



S. = ^ I d'xCcsix) (4.1) 



where /i is the regularisation scale that keeps the action dimensionless in = 3 — 2e. The path- 
ordered exponential is invariant under dilatations and the Lagrangian is covariant with weight 
Ac = 3, whereas the measure d'^x does not match this weight for d = 3 — 2e. This results 
in a non-vanishing variation of the action with respect to dilatations and special conformal 
transformations. The conformal Ward identities can be derived by acting on both sides^ of (4.1) 
with the generators of conformal transformations, see [46-48J. This leads to the Ward identities 

2?V f 

» (Wn) = -— / d''x{C{x)Wn) , (4.2) 

K'^iWn) = / d''xx-'{C{x)Wr.) , (4.3) 

for dilatations and special conformal transformations. Here the operators on the left hand sides 
act in the canonical way on the coordinates of the cusp points, 

D = 5^(a;. ■ 90 , ^'' = Y1 (2<(^^ " " ^'^0 • (4-4) 

i i 

We emphasise that thanks to the factor of e on the r.h.s. of fl4.2p and (4.3) it is sufficient to know 
the divergent part of the integrals appearing on the r.h.s. of those equations in order to obtain 
information about the finite part of (Wn)- 

The dimensionally regularised Wilson loop (Wn) is a dimensionless scalar function of the cusp 
points x^, which appear paired with the regularisation scale as x^jfj,"^. As a consequence (Wn) 



"^Note that the left-hand side of (4.1) is a function of the cusp points, whereas its right-hand side contains ah 
fields of the Lagrangian. 
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satisfies 

J2{x,-d,)-fi—j {Wr,) = 0. (4.5) 
This provides a consistency condition for the right hand side of (4.2). 

4.1 One- loop insertions 

At order N/ k we have a contribution from the contraction of the kinetic part of the Lagrangian 
insertion with the second order expansion of the Wilson loop operator, shown in fig. 5, 

(£(x)W^„)(i) = (£kin(x)H^„)(^) = ^ [dz^fe^^''{TT{A^dpA,){x)Ti{A^A,))('^ . (4.6) 

J Zi>Zj 

The direct calculation of the right hand sides of (4.2) and (4.3), yields a vanishing result as 



Figure 5: Lagrangian insertion contributing to the Ward identities at one loop, 
e — 7- 0. Thus we have 

D(iy„)« = 0(e) , and K'^(iy„)W = 0(e) , (4.7) 

in other words the conformal symmetry is unbroken for e = 0. As a consequence, the expecta- 
tion value of the Wilson loop is constrained to be a function of conformally invariant variables. 
Starting from the Lorentz invariants xfj the most general conformal invariants are the cross-ratios 

^ijkl ■■- -2-2- • 

In our case where neighbouring points are light-like separated, xf^^^ = 0, non- vanishing cross- 
ratios can only be written down starting from n = 6. The special conformal Ward identities (4.7) 
then imply that {Wn}^^^ is given by a function of conformal cross-ratios, 

(Wn)^'^ = gn {Ui.ki) , (W,)^'^ = const. (4.9) 

Since there are no non- vanishing conformal cross-ratios at four points, (^4)''^^ must be a constant. 

Let us now compare against the results of our one- loop computation of section 2, There, the 
constant on the r.h.s. of the second equation in (4.9) was found to be zero for the tetragon. 
Moreover, analytical investigations of certain symmetric contours and numerical investigations 
for non-symmetric contours show that geiuijki) is zero for the hexagon. As mentioned before, we 
expect that the result remains true for higher polygons, i.e. gn = 0. 
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(g) (h) (i) 

Figure 6: Two- loop contributions to the anomalous conformal Ward identity. The grey point in 
the pictures denotes the point x of the Lagrangian insertion in {C{x)Wn)- The diagrams of the 
second and third line cancel pairwise. 



4.2 Two-loop insertions 

At two loops there are several diagrams that contribute to the insertion of the Lagrangian into 
the Wilson loop, {C{x)Wn), that correspond to the kinetic term, the gauge field vertex, the ghost 
kinetic term and the ghost vertex in C{x). Those diagrams are shown in Figure 6, We do not 
display diagrams that vanish for kinematical reasons as at one-loop level. 

Just as at one-loop level, only diagrams giving rise to divergent integrals will contribute to 
the anomalous Ward identities. 

4.2.1 Insertion into the ladder diagram 

Let us consider the insertion of the kinetic term of the action into the ladder diagram as shown in 
Figure 6a, For the dilatation Ward identity these are exactly the two-gluon diagrams calculated 
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above, which are finite. For the special conformal Ward identity the integration is shghtly more 
comphcated, but the finiteness is easy to check for all contributions. Thus, this diagram does 
not contribute to the anomalous Ward identities. 



4.2.2 Insertion of the interaction term 

Next, we can contract the cubic order expansion of the Wilson loop with the vertex term of the 
Lagrangian insertion, as shown in Figure 6b, 

{C{x)W,)^^^ = {CU^)^{^f j rf'<;:nr {A.A^A,)) . (4.10) 
For the dilatation Ward identity, we trivially have 



/ 



2 



d'x{C{x)W,)'§^ = ]{W,)^2te. = '(j) ^ + Oie') , (4.11) 



e 



which is just the vertex diagram that was calculated in (3.6), up to a factor of i. 

The contribution to the special conformal Ward identity is more complicated. We have 

with 

Jo J \X\'^\X - Zi2\'^\x - Z32r 

= ^^^^{x2 + x,r + 0{e'), (4.14) 
e 

where the coefficient ln(2) was computed numerically to 10 relevant digits. The reason the pole 
arises was discussed in section 3.2. 

Details of this calculation can be found in appendix IB. 3. II Summing up all four contributions 
^321 > ^421, ^432, ^431 ^c arrive at 

j ^."'{rc.jwyg = 1 (4.15) 

4.2.3 Insertion of the kinetic term into the vertex diagram 

Furthermore, we can contract one gauge field of the kinetic term of the insertion with the Wilson 
loop and the other one with the 3-gauge-field vertex, leading to a diagram of the type displayed 
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in Figure l6cl, 



{C{x)W,)f]^ = (/:kin(x)^ fdz^;;X Tr (A.A^A,) (^z J d'^wC,^^ ) (4.16) 



A; \^ 1 



^a/.7g^-(Tr (A^dpA^) (x)Tr (A^A^A^) Tr (A5A.A,) («;)) 

Let us Wick-contract the kinetic term with A^{zj) (the two other contractions are discussed 
below.) We obtain 



d 

where G^y = t^yp{x — uY /{—{x — y) and where 



luaix - Zi,x -w) = e"^'^ 



Gay{x - Zj)df^G^fj{x - W) + Gaa{x - w)df^G^u{x - Zj) . (4.17) 



The only dependence on the insertion point x is in ly^j. For the dilatation Ward identity the 
integral / d'^xl^^ can easily be computed (for details see Appendix B.3.2) and effectively gives a 
propagator such that we have 

j d^x (£(x)1^4)S5 = 3^(1^4)™''*^^ = -32 + 0(e°) , (4.18) 

where a factor of 3 was included since the insertion can be in any of the 3 propagators of the 
vertex diagram and thus we get three times the same contribution. 

The contribution to the special conformal Ward identity is more complicated, since the inte- 
gration J d'^xx^Iijfj does not just yield a propagator. Performing the calculation, we find 

j d'^xx^{C{x)W,)^^ = (f )'^ln(2) l^xr + 0(6°), (4.19) 
where the coefficient ln(2) was computed numerically (details can be found in Appendix B.3.2). 



4.2.4 Insertions with gauge field and ghost loops 

The gauge-field-ghost-insertions in Figure 6 cancel pairwise. For the dilatation Ward identity 
the insertions of the three gauge-field and the gauge-field-ghost vertices as shown in diagrams 
6d, 6g are identical to the gauge field and ghost loop diagrams (3.12) and thus cancel. It is not 
necessary to perform the integration over the insertion point to see how the cancellation occurs 
and thus the contributions to the special conformal Ward identity cancel as well (for details on 
the cancellation see B.2). 
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Insertions of the kinetic term into the gauge field propagator as shown in diagrams 6e, 6h 
cancel as well, since the insertion is the same for both diagrams and thus the algebraic relations 
responsible for the cancellation in (3.12) remain unchanged. 

Inserting gauge field respectively ghost kinetic terms into the propagators inside the loop in 
diagrams 6f, 6i produce slightly more complicated expressions. Nevertheless they cancel as well 
as can be seen in a straightforward calculation. For the dilatation Ward identity the cancellation 
can be seen in an even simpler way by noticing that the integration over the insertion point x 
effectively yields a gauge field respectively ghost propagator. Thus the diagrams are identical to 
the ones in (3.12) and cancel. 



4.3 Anomalous Ward identities and generalisation to higher polygons 

Summing up the divergent contributions of (4.11) (4.18), and inserting them into the dilatation 
Ward identity (4.2), we obtain 

B{W,)('^ = - (f ) ' ln(2) 1 j + 0{e) , (4.20) 

where we have written the factor 4 as (X]^=i 1) ^'^ emphasise its origin from the sum of four 
vertex-type diagrams. Note that only the divergent part of the vertex- diagram was required here. 
Summing up (4.15) and (4.19), and inserting them into the special conformal Ward identity (4.3), 
we obtain 

//V\ 2 4 

K'^(H^4)(2) = _2 / -^j ln(2) ^xr + 0(e). (4.21) 
^ ^ i=i 

Let us now explain how these equations can be generalised from n = 4 to arbitrary n. In our 
two-loop computation, we found that the only diagrams contributing to (4.20) and (4.21) are 
those producing poles in e. The mechanism for how these poles are generated was described in 
section 3.2, see in particular Figure 4. It is clear that for > 4 cusps, the same type of vertex 
diagram will produce the divergent terms. Although those diagrams will depend on one further 
kinematical variable w.r.t. the four-point case, this dependence cannot change the (leading) UV 
pole of the diagrams. Since there are n diagrams of this type at n points, we expect 

B{W^)('^ = - (f ) ' ln(2) 1 j + (9(6) , (4.22) 

and 

K-(|y„)(2) = _2 _ ln(2) 5^xr + 0(e),. (4.23) 
^ ^ i=i 

We will now proceed to discuss the solution of these Ward identities and compare them to the 
result of the two-loop computation of the tetragon Wilson loop in section 3.4, 
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4.4 Solution to the anomalous conformal Ward identities 

Using D {^Ij) = '^x'ij it is clear that the most general solution to the dilatation Ward identity 



(4,2.2) is 



(2) 



ln(2) 



E 



+ fn 



X7 



X 



kl 



0(e) 



(4.24) 



Of 



where / is an arbitrary function of dimensionless variables and we recall that /i^ 
course, this is exactly what we expect from (4.5). 

The result for the special conformal Ward identity is more interesting. Plugging (4.24) into 
the special conformal Ward identity (4.21) and using K*^ ln(a;^;) = 2{xk+xiy , it is easy to see that 
the function /„ is allowed to depend on conformally invariant cross-ratios only, i.e. fn{x]j/xli) 
= giuabcd)- Therefore we finally have 



(2) = _ 



ln(2) 



E 



-Xi 



+ gn{Uabcd) + 0{e) 



(4.25) 



In the four-point case, there are no non- vanishing cross-ratios, and therefore in particular must 
be a constant. This is in agreement with (3.13) and thus represents an independent check of the 
direct perturbative computation, including its finite part (recall that deriving the Ward identity 
does not rely on the finite parts of the direct perturbative computation). So, even though the 
result for the vertex diagram (3.11) was obtained numerically, its functional form is an analytical 
result, since we know the analytical expression for the ladder diagram and the sum of vertex and 
ladder diagram through the solution to the anomalous conformal Ward identity. 

5 Two loops: ABJM theory 

Here we explain how the results are modified in ABJM theory. We use the Wilson loop operator 
proposed in [6J 



{W{A,A)) 



2N 



TrPexp i (b A^dzi" + TrPexp i (p A^,dz^ 



(5.1) 



Note that the sign(s) in the exponent(s) in (5.1) are correlated to corresponding signs in the 
Lagrangian by the requirement of gauge invariance, see Appendix A. 



5.1 Gauge field contributions 

In ABJM theory there is a second copy of the gauge field with opposite sign in the Lagrangian 
(A. 20). Up to a sign, the gauge field contributions for both gauge groups are identical at one 
loop. 



(5.2) 
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due to the different sign of the propagator for the second gauge field, {Af^A^) = — (A^A,^). Thus, 
at one loop the diagrams cancel. This does not differ from the result of pure Chern-Simons 
theory, since the expectation value at one loop vanishes, as we found in section 2 for n = 4, 6. 

At two loops, however, the sign has no effect, since the two-gluon diagram contains an even 
number of propagators and in the vertex diagram we have to take into account the sign of the 
interaction term as well. Therefore, the two-loop diagrams are identical 

{W)^I^ = {W)f (5.3) 

Thus, up to two loops, the expectation value for pure gauge field contributions is the same in 
ABJM theory and Chern-Simons theory 

{W{A, A))„ fields = {W{A))cs . (5.4) 
5.2 Matter contributions 

In pure Chern-Simons theory the one-loop correction to the gauge field propagator is zero, since 
the contributions of gauge fields and ghosts exactly cancel against each other, see (3.12). 

In ABJM theory we have to take into account fermionic and bosonic matter in the loop. This 
gauge field self energy has been calculated in [49,50,6J and the corrected propagator readsi, 



'"'^ ^ \k J 8 r(d-l)7r'^ I r(2-f) (-a;2)rf-2 "^"M ^(3 - f ) 4 (-x2)rf-3 



(5.5) 



for details see Appendix C. We can drop the derivative term in (5.5) (it would not contribute to 
the gauge-invariant Wilson loop) and instead use the propagator 

Gi;>W^-^(f)^■-^(f-l)^^, ,5.6) 

which up to two small differences is the tree level Af = 4 SYM gluon propagator. The first 
difference is a trivial prefactor, and the second is that since we are at two loops, the power of 
is 1 — 2e here, as opposed to 1 — e in the one-loop computation in A/" = 4 SYM. Thus 
it is clear that the results will be very similar to the expectation value of the Wilson loop in 
A/" = 4 SYM. The corresponding calculation in A/" = 4 SYM was carried out in [13j and we briefly 
summarise the results. 

As in A/" = 4 SYM we have three classes of diagrams shown in figure 7. Diagram 7a van- 
ishes due to the light-likeness of the edges, whereas 7b yields a divergent, and 7c yields a finite 



^Recall that we absorbed the regularisation scale into the couphng constant: fc — > /i ^^/c 
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(a) 



(b) 



(c) 



Figure 7: Examples for the three classes of diagrams involving the gauge field self energy: di- 
agrams with (a) a propagator connecting the same edge (vanishing), (b) propagator stretching 
between adjacent edges (divergent), (c) propagator stretching between non-adjacent edges (fi- 
nite). These diagrams have the same structure as the 1-loop diagrams in A/" = 4 SYM. 



contribution. We have 



(2) 



matter 



(5.7) 



Zi ^ Zj 



i>3 
2 



(jl ((4-™^" 



\2\2-d 



where lij = J dsidsjPi ■ pj{ — {zi — 

There are four diagrams of the type shown in fig. 7b, It is sufficient to compute one of 
them, as the others can be obtained by the replacement i ^ i + 1. Setting e.g. i = 2, j = 1, we 
have 



1 



hi = -2(^^13 



2 \3-Q 



ds2dsi- 



1 



-X 



2 \2e 
13 J 



[S1S2) 



vl-2e 



Furthermore, there are two finite diagrams Ii+2,i of the type shown in fig. 7c. Setting (i = 3 and 
taking i = 3, j = 1 we have 



'31 



- / ds-idsi 



^2 I ^2 
■^13 "T "''24 



Taking the sum over all contributions we obtain 



InM ^ 1 + vr^ 

X24 



i>j 



2 \2e 



2 \2e 
24 J 



+ 21n^ 



X 



13 



X 



24 



+ 2tt' 



(5.9) 



(5.10) 



and thus the full matter part is 



(W4) 



(2) 



matter 



1 fN 

"4 It 



-xj^ine^'^fi 



2 \2e 



-x|447re]V]^_2ij,2 



X 



13 



24 



-7i^ + 0{e) 



(5.11) 
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where we have restored the regularisation scale /i^. The full result is obtained by adding the CS 
part (3.13) to (5.11) and the result can be rewritten in a form in which the terms cancel 



(2) 

ABJM 
,/2 



2 r 



'2^2 \2e 



+ 



21n^ 



13 



24 



- const. + 0{e) 



(5.12) 



where /i' = Svre'^^/i^ and const. = 81n(2) + 20 In (2) + 27?^ — og and where ag is given in (B.27rt 
and fits the value oe = -^ir^ + 161n(2) + 81n^(2) such that const. = fyr^ + 121n^(2) - 81n(2).l 
This is the result quoted in the introduction (1.6). 
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A Conventions 

An n-sided polygon can be defined by n points Xi {i = 1, ...,n), with the edge i being the line 
connecting Xj and Xj+i. Defining 

Pi = <+i - < (A.l) 
and parametrising the position on edge i with the parameter Si G [0, 1] we have 

^f(.,)=xr + pf.,. (A.2) 

Furthermore, we use the notation 

e{p,q,r) = e^upP^q^^r^ and Si = I - Si. (A. 3) 

One can easily check that with the definition x^- = xf — 

2 Xij ■ Xfnn Xj^^ + Xj^ ^im -^jn (-^-4) 

We consider 3-dimensional Minkowski space with metric 

T]^, = diag (1, -1, -1) . (A.5) 
Using (A. 4) we can rewrite the scalar products 

2 2 2 2 / \ \ 

Furthermore, using the definition (A.2), one can easily show that 

(^i "~ Zj) = X^jSiSj + X^_^_ijSiSj + Xj^j^iSiSj + X^_^ij_^^iSiSj . {■^■'^) 



^Note that the log(2) term has a different transcendentahty than the other terms. It stems from the vertex 
diagrams of the pure CS part (3.13) of the Wilson loop. 
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A.l C herns- Simons theory and the Wilson loop operator 

The Chern-Simons action 



Scs = ^j d^x e^""" Tr f^^^d^A, - ^A^A.a}^ (A.8) 



(A.9) 



and the Wilson loop operator 

W\C\ = ^TrPexp fyi j A^dz^ . 

are invariant'^ under SU{N) gauge transformations 

A^^A'^ = g{x) {A^ + zd^) g-\x) , g{x) G SU{N) (A.IO) 

if the path ordering^ is defined as 

V{A,{z{s))A,{z{s'))) = A,{z{s))A,{z{s')) for s > s' , (A.ll) 
r{A,{z{s))A,{z{s'))) = AAz{s'))A^{z{s)) for s < s' 

and s G (0, 1) parametrises the path along the curve C. The path ordered exponential in the 
Wilson loop operator then has the expansion 

Pexp (^i A^dz''^ = Inxn + i dz>;A^ + {if jf rfzf j ' dz-^A^A^ (A.12) 

+ {tf j( dz^ J ' dz^ J ' dz'^A^A.Ap 

+ {{)^jdz'^J dz'^ j dzlj dz1A^A,ApA^+ .... 

Quantising the theory with the standard Fadeev-Popov procedure yields the gauge fixing and 
ghost action 

5g., = ^jd''x Tr Q {d^A^f + c [d^-D,) ^ (A.13) 
where D^c = (9^c + c]. In Landau gauge = 0) the gauge field propagator reads 



p 



{{A^) (x) (A,),, (y)) = 6u6,kT ^^up— (A.14) 



I 2 



^More precisely, the action is invariant under infinitesimal transformations 5(2;) = 1 + ia(x) and transforms 
like 5* ^ S' + (2Trk)SS, where SS = J d'^xe'^'PTr {{df,g-^)g{d^g-^)gidpg-^)g) under finite transformations. 

Since dS takes integer values, exp(i5') is invariant under large gauge transformations for fc G N. 

*There are different conventions on the path ordering. They are equivalent, it is however important to choose 
the sign in the Wilson loop such that it is gauge invariant. Reversing the sign in the exponent of the Wilson loop 
does not just reverse the integration contour, due to the path ordering. 
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where we have rescaled the couphng constant k /j, "^^k and restore the dependence on the 
regularisation scale /x only in the final results. The ghost propagator is 

(c(x)c(y)) = 5u5,,- I > I — . (A.15) 



TT 2 



Note that the gauge field propagator is related to the ghost propagator by 

((A,)^^. (x) {A,),, (y)) = ^e,,^dP{c{x)c{y)) , (A.16) 

which can be used to see the cancellation of gauge field and ghost loop contributions to the 
one-loop gauge field propagator in a simple way. 

The different conventions found in the literature on Chern-Simons theory deserve a short 
comment. One can show, that 



k 



S = — I d'^xTi I A^d.Ap - -sA^A.Ap ] e^'^" , (A.17) 



W{C) = ^TiVe^^i^s j A^.dz"^ (A.18) 

are invariant (In the sense mentioned above) under gauge transformations 

A^^A'^ = g{x) (^A, - -d)j g-\x) (A.19) 

where s is a real or imaginary parameter. I.e. the sign in the Wilson loop and the Lagrangian are 
correlated through gauge invariance. Taking a hermitian gauge field (^4^)"'' = we can choose 
s = i, which corresponds to the choice used throughout this document. 



All other conventions found in the literature can be obtained by rescaling the gauge field 
A^ — )• —A^, A^ — )■ iA^ etc. Note however, that this changes factors in the Lagrangian, the gauge 
transformation, the covariant derivative and the Wilson loop. A sign difference in the Wilson 
loop only may also be due to a different definition of the path ordering (A. 11). 

A. 2 Lagrangian of ABJM theory 

The action of ABJM theory is 

i^ABIM = <ScS + <Sg.f. + Scs + <Sg.f. + ^matter (A. 20) 
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where S is obtained from S by replacing with the gauge field in the anti-fundamental repre- 
sentation Afj_ and letting k — )■ —k. Explicitly, we have 



2 . 

Tr^A^d^Ap — —i Ap_A^Ap) — Tr^A^d^Ap — —i A^A^Ap 



2 

—I 

3 



S'gf + Sgf 



k 

4:71 



d X 



(Tx 



Tr(D^ Ci D^C^) + z Tr(7/'^ p ijj) 



(A.21) 
(A.22) 
(A.23) 



The field content consists of two U{N) gauge fields {Ap)ij and {Ap)^, the complex fields {Ci)--- 
and {C^)ij^ as well as the fermions {ipi)'^ and (V^^)j- in the (N, N) and (N, N) of U{N) respectively, 
I = 1, 2, 3, 4 is the SU{4)fi index. We employ the covariant gauge fixing function dpA'^ for both 
gauge fields and have two sets of ghosts (c, c) and (c, c). S'int are the sextic scalar potential and 
i[r2(j2 Yukawa type potentials spelled out in [5]. The covariant derivative is given by DpC = 
dpC + i[Ap, .]. It's action on , is given in [6J and not needed here. 

B Details of the two-loop calculation 

B.l Vertex diagram 

Using e{pi,p2, Z12) = €{p3,p2, Z32) = the first line in (3.8) can be rewritten as 

{d-2)-' 



/32I = / Si,2,3e(pi,J92,5zi)e(P3,P2,5x;3) dw 



\w\'^\w — Zi2\'^~'^\w — Z-i2Y~'^ 



(B.i: 



We begin by introducing Feynman parameters 

1 



wY'\w — Z\2Y'~'^\w — -232!'^"^ 



where 



and 



(-(«;- /3lZl2-/33^32)^ + /^)^''"'^/' 

d/3ici/32c//33(/3i/32/33)^'~'^/'-'/925(5^ A - 1 



r(f -2) 



r(f)r(f-ip 



(B.2) 
(B.3) 



A = 2f3Mz,2 ■ ^32) - - 4/^3^3 • (B.4) 

Shifting the integration contour w ^ I = w — f3iZi2 — fizZz2 we have a standard integral 



dH 



[/2 - AY 



r(n f I 



r(n) \A 



Thus we get 



Cl 



'321 



id -2) 



id Si,2,3e(pi,P2,<9^i)e(P3,P2,<9^3)^ 
25 



-2 



(B.5) 



(B.6) 



where ci = m^r [d — 2)/ F (3(i/2 — 2). Evaluating the action of the derivatives and abbreviating 



a;^3 = s, a;24 = t we obtain 



/321 = c^st J^' dh,^,,sd'/3i,2AP,/3,/3sY'-'^/'6{J2 A - 1) - 2^^^ms,ss{s + t)) 

(B.7) 

where C2 = irr'^^'^F^d — 1)/ {8F^(d/2)) and both terms are separately symmetric under s t. 
Performing the change of variables 

(3i = xy, (32 = xy, f^s = y , XI ^* = ^ x,yE [0,1] (B.8) 

i 

with Jacobian y we can rewrite (B.38) in a form where all integrations range from to 1, 



Ia = C2st I ' dhdxdy ^'^'^l" , (B.9) 

y 

I -\d _ d-2 _ 

[xy)2x 2 S1S3 



Ib = —2{d — l)c2st{s + t) d^sdxdy 

Jo 



where 

Ay = - {sxsi{s3y + S2xy) + tysz{xsi + S2x)) , (B.IO) 
and /321 = /a + /b- The integral I a is divergent as si, S3 — )■ 0, see also fig. 4] 

B.1.1 Numerical evaluation using the Mellin-Barnes method 

In this section, we switch from the Feynman parametrization in equation (B.9) to a Mellin- 
Barnes representation, as the latter is very convenient to perform a systematic expansion in e. 
An introduction to the Mellin Barnes technique can be found in [51j. 

In the first step the sum in the denominator is transformed into an integral over a product 
of terms. Since the denominator in (B.38) consists of a sum of four terms, we will introduce 3 
Mellin parameters zi, Z2,Z3. By repeated use of the Mellin-Barnes representation 

1 11 rP+ioo Y 



{X + Y)^ F{X) 2m 
where — Re(A) < (3 < 0, one obtains 
{2mfF{X) 



/ -^^F{z + \)F{-z)dz, (B.ll) 



[a + h + c + dy 



j dzi,2,3a''F'c''d-^-''-''-''F{-zi)F{-Z2)F{-Z3)F{X + zi + Z2 + Z3) (B.12) 
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where the real parts /3j of the integration contour have to be chosen such that the arguments in 
all r functions have positive real part. Applying (B.12) to the denominator of I a (B.9), we can 
rewrite I a as 

T ■ - ^ j d~Z,,2,3n-Zl)n-Z2)n-Zs)r{d -l + z, + Z2 + ^3)(-s)^i+-^ + l(-t)-'^--i-^+2 

(B.13) 



r{d 



I 



where dz = (27rz) ^dz. The integrals over si, S2, S3, x, y can be carried out using 

s'^-\i - Sif-^t = r{a)r{b)/r{a + b) , (b.m) 



and we arrive at 

C2 



J dh,2A--^r^''^\-tr'~'''~''^'n-^i)n-^2)n-zs) 



(B.15) 



^ r{d 

X r{zi + i)r{z2 + i)r(z3 + i)r{d -i + zi + z2 + z^)r {-d/2 -Z2 + 1) r{zi + ^2 + 1) 
X r{-d -Z1-Z2 + 2)r {-d/2 -zi-zs + i) r{-d -zi-z2-z^ + 2)r {d/2 + Z1 + Z2 + z^) 
X [r (2 - d/2) r{-d -Z2 + 3)r{-d -zi-z^ + 3)r{zi + Z2 + Z3 + 2)]-^ . 

Recall that we investigate the kinematical region where s,t < 0. 

One can see that this integral is divergent as e — by noticing that for e = it is impossible 
to choose the integration contours such that all poles of r{... + zi) are to the left of the integration 
contour and all poles of r{... — zi) are to the right of the contour. ^ The reason is that the 
poles of r{zi + ^2 + 1) and r{—d — zi — Z2 — z^ + 2) = r{—zi — Z2 — I + 2e) "glue together" at 
zi = —Z2 — 1 for e = 0. However, one can find allowed contours for e ^ 0. 

By shifting the contour left to the pole at zi = —Z2 — 1 we pick up a residue. The factor 
of r{—zi — Z2 — 1 + 2e) evaluated at the residue, results in a divergent factor of r{2e). The 
remaining integral over the shifted contour yields a finite contribution. 

The steps of shifting contours and taking residues have been automatised in [52j and we used 
this package to systematically extract the pole terms. Applying this procedure to (B.15) and 
expanding in e yields 3 integrals: 

1^1^ = C2 j ^(i + 21og(-s)-^i(^0)/i(^i) + C?(e) (B.16) 

lf=C2 j ^0Q + 21og(-s)-^?2(^l,^3))/2(^l,^3) + O(6) (B.17) 

lf = cJ'-^^{T^''''h{z.^.z^ (B.18) 



{2mf \t 



^This is necessary in order for the previous steps to be well-defined. 
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We do not specify the values of the real parts /3j as well as the functions fi, gi here, which are 
lengthy expressions of products of F functions and can be obtained automatically by expanding 
(B.15) with the help of [52j. Adding up the divergent part of (B.16) and (B.17) we get 



where by numerical integration one finds 

«i= / / ^^/2(^i,^3) = 8.710344 ±10-''^47rln(2) = 8.710344361... (B.20) 

J 27r2 J [Z-Kiy 

accurately approximated by our analytic guess. Further numerical evaluation of the finite part 
oiI^^\l^2\lf yields 

/yi,fimte = ^2 {ai In(-s) + Oi In(-t) + a2 \y? + . (B.21) 
where ai has the same value as above and 

02 = -0.84 ±0.01 04 = -14.375216465 ± 10"^ (B.22) 
Numerical analysis for Ib suggests 

/b = C2 (aaln' (^) +05) (B.23) 

where 

as = 3.97 ±0.01, 05 = 40.620843911 ± 10"^ (B.24) 
Adding up 02, 03 we obtain 

a2 + a3 = 3.136 ±0.02 ^ vr. (B.25) 

Thus suggesting 

Ia + Ib = C2 ± ai (In(-s) ± In(-t)) ± (02 ± 03) In (^) ± (04 ± 05) ± C(e)) (B.26) 
~ C27r l^iM^ + 4 in(2) (In(-s) ± In(-t)) ± In^ (^-^ +a^ + 0{e)^ 

where 

ag = (04 ± a5)/7r = 8.354242685 ± 2 ■ 10"^ . (B.27) 

The constant fits the value ^ -fvi^ ± 161n(2) ± 81n^(2) = 8.35424273... . The result can be 
rewritten in the form 

Ia + Ib^ C27r (^21n(2)-^^^^£l^-^^^ ± In^ ± ag ± 0(e)^ (B.28) 
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B.2 Gauge field and ghost loops 



It is well known [35j, that the contributions of ghost and gauge field loops to the gauge field self 
energy cancel. We briefiy review the cancellation of the gauge field and ghost loop corrections in 
the Wilson loop, since from this it is easy to see, how the cancellation for the insertions in the 
conformal Ward identities takes place. 



B.2.1 Gauge field loop 

The gauge field loop-diagram arises at second order in perturbation theory 

^^^^giuon-ioop ^ ^^Tr (- f dz^dz^A.A,) (-i) (A I d'xe'^^^Tr (^^A.A.A,)) ') (B.29) 

^ ^dz^dz-^ j dVye^^'<e'''^i:i{A^A,)i:i{A^ApA,){x)i:i{AsA,A,){y)) 



= cs 
where 



1 1 / ^ ^ 



Taking into account that A^, A^, give 3 identical contractions with one of the vertex terms and 
that we can contract them either with the x- or ?/-dependent vertex, we get a symmetry factor 
of 3 • 3 ■ 2. The remaining contractions of the gauge fields are dictated by taking into account 
only planar diagrams. Thus we get 

(iy„)s"°P = Sc^jdztdz^; j d''xd''ye''P'^e'''^{A,A^){A,As){ApAr){A^A,) . (B.31) 

To proceed, we recall the relation (A. 16) between gauge field and ghost propagator and write 

e"^^e^--(A^A.)(A,A.) = ie"^^e^'^-e^.«e,.,9,^(c(x)c(y))9,«(c(x)c(y)) 

= \d'M^)c{yWy{c{x)-c{y)) (B.32) 

where we used e'^^'^ e^"'^ e^rK^'yup = — {UkVI + V^V^k) ^"^^ d^F^x — y) = —dyF{x — y) in the last 
step. 



B.2. 2 Ghost loop 

The ghost loop diagram arises from contraction of the second order perturbation theory expansion 
of the gauge-field-ghost vertex term 

^^^^ghost loop ^ i_(Tr (-^ rfzfrf^^A.A.) (-1^ (A I d'xTTidPczlA^c])^) (B.33) 
= 2c8 (fdz^dz-^ I d''xd'^y{Ti{A^A,)i:i{dPcApc)Ti{dlcA„c)) 
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where cs ist the same as defined above and the factor of 2 is due to the fact that the evaluation 
of the first fine yields two identical planar diagrams that are kept and two identical non-planar 
diagrams that we drop. Contracting either with the x- or y-dependent vertex, we get a 
symmetry factor of 2. There is only one way for the remaining contractions and thus we get 

V ghost loop __ /I i . U , 1. / Jd^,d 



loop ^ ^ ^^,^^u J d'':,dMAM{AuAp)dPMynx))d'^{c{x)c{y)) (B.34) 

where a factor of —1 due to the anti-commuting ghost fields in the loop was taken into account. 
Summing up (B.31) and (B.34) we get 

(PVn)™ field loop ^ ^ly^^ghost loop ^ g _ j^g 35^, 

The same relation (A. 16) can be used to show the vanishing for the dilatation and special con- 
formal Ward identities. 



B.3 Conformal Ward identity 
B.3.1 Insertion of the interaction term 

We can rewrite (4.13) as 



{d 



Introducing Feynman parameters, changing the integration variable to / = x — (3iZi2 — /33Z32, using 
the same notation as in app. B.l, integrating over / and evaluating the action of the derivatives 
yields 

/321 = J d's,,2,3 J d[(3],(^{P,z, + P2Z2 + /33Zsreip,,p2,dMP3,P2,dJ^ (B.37) 

{2-dy 



+ e(pi,p2,P3)2/3i/33e„/p^ iPiPiSi + P3P3S3) 



The last term can be shown to be finite and the first term is very similar to the vertex diagram. 
Evaluation of the derivatives as in B.l yields 

^321 = C2st f rf'si,2,3Ci'/3l,2,3(/3l/32/33)^'~'^/'5( V A - 1) (B.38) 
Jo i 

(/3i2;i + P2Z2 + P3Z3y ( ^ - 2^^^^/3i/33SiS3(s + t) j + finite . 
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It can be shown, e.g. using the Melhn Barnes technique as in B.l, that all divergent contributions 
are due to the first term. We have the following divergent contributions: 

St J dhd[f3W2(33 = i a4 + O(e0) , (B.39) 

St [ d'sd[PW2P3 (^) =-bixl^ + x^) + Oie') , 



St J dhdlPW^Ps = ^ (a 4) + 0{e') . 

Numerical evaluation yields 

a = 1.8562 ±0.0001, 6 = 2.4989 ± 0.0001 . (B.40) 

To good accuracy we find 

a + b = 4.35517 ± 0.0002 ^ 27rln(2) = 4.35517... . (B.41) 
Summarising I'^2i then reads 

4i = + b){x2 + xsY + 0(e°) ^ ^^^^^{X2 + xsY + 0(e°) . (B.42) 

B.3.2 Insertion of the kinetic term into the vertex diagram 

For the kinetic insertion into the vertex diagram we have (4.16) 

(^(^)W^4)SeJ = (f ) ' ^ (^) / j rf<;S^''^^^- - ^)GA^k - w) (B.43) 



+ cyclic(/i, i/,p;Zi,Zj,2;fc) 



where G^,y{x - y) = e^^p— ^^^^ and 



Iu(t{x — Zj,X — w) = e"^^ (Cauix — Zj)d'^p^G^ia{x — w) + Gaa{x — w)d^^^G^iu{x — Zj))^ . (B.44) 

and the two other contractions are contained in cyclic(/i, p; Zi, zj, Zk). For the dilatation Ward 
identity the integration over x can be performed by introducing two Feynman parameters. The 
result simply yields a propagator 

d Xli,^ = ——Jjr^uaifi- - —J = ~ 7-. (d\ Gvo-{Zj — w) , (B.45) 
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In the case of the special conformal Ward identity the integration is a httle more involved. The 
integral over d'^x can be solved by introducing Feynman parameters and after some algebra one 
finds 

j d''xx\C{x)W^)fl = 2c^c^ j dz'^;;;^ j rfV^^e^^^e^.^e^,,^ (B.46) 

1 \ ( {zj+wY 



{2r,-% + 4d-) - d-d, 

V \{{zj-wY)2 y^^Zj-wy)-2 

GpsG^,r + cychc(/i, z/, p; Zi, Zj, Zk) , 

where all derivatives are taken with respect to zj and C4 is a constant obtained through integration 
over d'^x. Inserting the propagators we can write this in a form convenient to solve the integral 
over d'^w 

2C3C4 ^ dz^^';;fe'''-e''^^e^,^e^^^e,s^e^redfdld^ {{'2v''%,(s + 4dJ)Jj - djd^jj}) + cychc (B.47) 
where the integrals read 



J = / d'^w 5 5 ^ (B.48) 



' {2-dyJ ((^. _^)2)f-2((^_y^)2)|-l((^^,_^)2)f-l 

The integrations over w can be performed by introducing three Feynman parameters and we 
get 



where 



/ dW,, = ^ dPdMM {^Pi - (B.50) 

and C5 is a constant obtained by integrating over w, the product C3C4C5 is explicitly given below. 
The expression for A is the same as in (B.4) For the cyclic permutations we can use the same 
expression, replacing the measure with (i[/3]3,i respectively (i[/3]3,A: , i-e. exchanging with 
respectively {3^^ in (B.50). 

All three contributions can then be written as 

2C3C4C5 j rf<;Je''^^e"''^e-ya^ [e^^i^e^sx^predldld^ {{2^''%^ + 4dt)Ji - dfdp.J^) (B.51) 
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and where 



2C3C4C5 = (^-J r(rf-4) (B.52) 

We can evaluate the derivatives and contractions with the computer and find that the non- 
vanishing contributions have the structure 

2 



(B.53) 



i>j>k ' 



where Iijk,p are lengthy terms proportional to l/A"^. 

For the conformal Ward identity we are only interested in the divergent part of the above 
quantities, which can be automatically extracted with the Mellin-Barnes technique. We find that 
all terms vanish except for i ^ j ^ k. Specialising to the case i = 3, j = 2, /c = 1 we find 



/3i2,-.-i = 0{S) , (B.54) 

/ /3V^+i = ^(6ix^ + &24) + 0(e°). 
Numerical evaluation of the integrals yields 

ai = 0.3465735 ± 10"^ ^ ^ ln(2) = 0.3465735... , (B.55) 
02 = 0.3465735 ± 8 ■ 10"^ ^ ^ ln(2) = 0.3465735... , 

= -0.8664339 ± 14 ■ 10"^ ^ ln(2) = -0.86643397... , 

62 = -0.8664339 ± 10^^ ^ ln(2) = -0.86643397... . 

Adding up the results, summing over all four diagrams and taking into account the corresponding 
prefactors we get 

C One loop gauge field propagator in ABJM theory 

Here we review the calculation of the one-loop correction to the gauge field propagator, see 
also [6j. We have fermionic and bosonic contributions in the loop and thus 

g^^Ap) = g^.'j'Kp) + Gil^'\p) (C.l) 
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where 



and 



d^k Tr (7^(^ + ^)7.^) 
{27ry k'^{p + kY 



We use the DRED scheme for Dirac matrix operations as well as for Levi-Civita tensor contrac- 
tions, i.e. we work in strictly d = 3 to obtain scalar integrands and only then continue the loop 
momenta to d-dimensional space to perform the integrals in d dimensions. This scheme has been 
shown to respect the Slavnov- Taylor identities up to two loop order in [35j. 
Then we have 

Tr (7^(]/ + V)7i.V) = 2 {-Ti^,u{p + k)-k + 2k^k^ + p^k^ + p^k^) . (C.4) 

The last two terms can be dropped, since they vanish when contracted with (C.2). The same is 
true for terms proportional to p^^Pu in the bosonic term. 
Summing up all remaining terms we get 

Introducing Feynman parameters, we have 

.r.I 2eo r f d'^^ k ■ {p + k) ,^ ^, 

+A'i;,-2,w I do: J (C.6) 
where A = —aap^. Then, we shift k = 1 — ap and drop terms linear in 

Using the standard integrals 

d'^l f _ i rfr(l-f) f d'^l 1 i r{2 



{CI 



and ex^^,exu5 = Vf^uVf^s - Vf^sVf^u we get 
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The standard formula.^ j 



(C.IO) 



leads to the Fourier transform of (C.9) 

'2ny N6jr{l-I)r{lf i^^^^f f rid -2) r,,, I rid -3)1 



rid-1) \ri2- l)i-x''Y'^ ^ \^ r(3 - f ) 4 (-x2)'^-3 
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